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ON  OPTIMUM  INSPECTION  SCHEDULES 


by 

S.  Patrick  Koh 


Abstract 

This  paper  treats  the  problem  of  determining 
the  minimax  inspection  schedule  for  detecting  failure 
of  a  component  or  system  when  inspections  have  a 
cost  and  cost  of  failure  is  proportional  to  the  length 
of  time  between  failure  and  detection.  The  mini- 
maxing  is  done  with  respect  to  all  failure  distribu¬ 
tions  having  a  given  mean. 


This  research  forms  part  of  a  doctoral  dissertation 
to  be  submitted  for  approval  by  the  Department  of 
Industrial  Engineering  and  Operations  Research, 
Columbia  University.  The  research  was  carried  out 
under  the  supervision  of  Professor  Cyrus  Derman. 
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ON  OPTIMUM  INSPECTION  SCHEDULES 


1.  Introduction 

Consider  a  certain  component  of  an  operating  system.  The 
system  must  operate  for  T  units  of  time.  T  may  be  finite  (the 
finite  horizon  case)  or  infinite  (the  infinite  horizon  case) . 

The  component  can  fail  at  a  random  time  Y.  However,  if  it 
fails  there  is  a  cost  incurred  that  is  proportional  to  the 
time  between  failure  and  its  detection.  Failure  can  only  be 
detected  if  the  component  is  inspected.  However,  there  is  a 
cost  for  each  inspection.  Our  interest  concentrates  on  the 
problem  of  determining  an  inspection  schedule  x  that  reconciles 
the  two  types  of  cost.  Whatever  occurs  subsequent  to  the 
detection  of  a  failure  or  to  the  completion  of  the  system's 
mission  at  time  T  is  not  considered. 

In  particular,  if  x  =  (Xq,x^,  . . .  ,x^^^)  ,  Xp=  0,  =  T 

in  the  finite  horizon  case,  or  x  =  (Xq  ,Xj^ ,  .  .  . )  ,  Xq=  0  in  the 
infinite,  horizon  case  is  a  given  schedule  of  inspection,  F 
is  the  distribution  function  of  Y,  c  is  the  cost  of  an  inspection, 
and  v  is  the  cost  per  unit  time  of  an  undetected  component 
failure,  then  the  expected  cost  associated  with  the  schedule  is 

C{x,F) 

"  ^r+1 

=  1  /  {(r+Dc  +  v(x  -  t)  }dF  +  (n+1)  c  (1-F  (T)  ) 

r-0  X  ^  ^ 

r 

,  for  T  <  “,  n  <  * 

=  Z  f  {(r+1)  c  +  v(x  .,-t) ydF(t) 

r=0  x^  ^ 

,  for  T  =  “>,  n  =  w. 


When  F  is  not  completely  known  the  problem  is  to  find  x 

to  minimize  U(x)  where  U(x)  =  sup  C(x,F),  and  ^  is  the  class 

Feff 

of  possible  distribution  functions.  Assuming  T  <  »,  Derman[2], 
originally,  considered  this  problem  for  the  class  of  all 
distribution  functions  for  non-negative  random  variables, 
allowing  for  the  possibility  that  detection  of  a  failure  at 
an  inspection  is  uncertain.  Explicit  formula  for  the  optimal 
schedules,  called  minimax  schedules,  were  obtained.  Roeloffs 
[5],  assuming  certain  detection,  obtained  minimax  schedules 
for  the  case  when  f.i  is  the  class  of  distributions  having  a 
known  quantile.  Kander  &  Raviv  [4]  assumed  Roeloffs 's  case 
with  the  added  assumption  that  all  distributions  in  ^  have 
increasing  failure  rates,  Beichelt  [1]  generalized  Derman's 
results  to  the  case  that  the  failure  costs  were,  rather  than 
being  linear,  an  increasing  function  of  the  time  between 
failure  and  detection.  Beichelt,  also,  considered,  in  the 
infinite  horizon  case,  with  increasing  failure  costs,  the 
minimax  schedule  when  Vf  =  consists  of  all  non-negative 
distributions  F  with  a  given  mean  p .  In  this  case,  he  proves 
that  an  optimal  schedule  has  equal  intervals. 

We  consider  for  the  finite  horizon  case  the  problem  of 

f  c 

determining  the  minimax  schedule  when  .  in  Derman  [2], 

the  observation  that  tf  can  be  reduced  to  the  class  of 

degenerate  distributions  (or  one  point  distributions) ,  is 
exploited.  Analagously,  we  employ  a  theorem  due  to  Hoeffding  (31 
that  enables  us  to  reduce  to  ^  ^  the  class  of  all  two-point 
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distributions  with  mean  u.  Results  are  obtained  that  permit 
the  numerical  determination  of  the  minimax  schedule.  We  also 
obtain  Beichelt's  infinite  horizon  schedules  by  the  limit  of 
minimax  schedules  as  T 


2.  The  minimax  schedule  where  only  the  mean  of  F  is  known. 

We  are  interested  in  the  problem  of  determining  the 
minimax  schedule  when  only  the  mean  of  F  is  known.  For  the 
infinite  horizon  case,  Beichelt  [1]  showed  that  the  minimax 
schedule  is  a  strictly  periodic  schedule  with  the  equal  interval 
space  6  ,  6  =  ' 

We  focus  on  the  problem  for  the  finite  horizon  case. 

We  start  with  this  problem  by  considering  the  supremum  expected 


cost  for  a  given  schedule  x,  that  is  sup  C(x,F).  Hoeffding  [3] 
proved  sup  C(x,F)  =  sup  C(x,F).  In  section  2.1,  we 


show 


0<i<m 


m+l< j  <n 


sup  {G. (x) ,  G. . (x) }=  sup  {C(x,F) } 

_•  1  11  -  O 


where  <  u,  >  u. 


G^(x)  =  (r+l)C  +  v(x^^j^-x^) 


,  r-0,1,2,  .  .  .  ,n,  Xq-0,  T 


X.  -  p 


p  -  X. 
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In  section  2.2.,  we  prove  a  key  theorem,  which  gives  a 
necessary  condition  of  the  minimax  schedule,  that  is 

(2.1)  -  G^(x)  _ 

’^r+l  “  ^r 

r=0 »1 ,2 , . . . ,n  for  some  Z  ^  0.  We  call  the  schedule  satisfying 
(2.1)  a  schedule  with  equal  average  increment  (SEAI) .  The 
supremum  expected  cost  of  a  SEAI  has  a  simple  form,  Gq (x)  +  . 

In  section  2.3.,  we  convert  this  problem  to  a  problem  of 
determining  the  minimal  point  of  a  continuous  function  with  a 
single  variable.  Under  some  certain  conditions,  the  function 
is  oiecewise  convex. 

In  section  2.4.,  we  show  that  the  minimax  schedule  for 
the  infinite  horizon  case  is  the  limit  of  the  minimax  schedule 
as  T  -*■  «>. 

In  section  2.5.,  when  T  <  an  algorithm  for  computing 
minimax  schedule  is  given. 


2.1.  Supremum  expected  cost  of  an  inspection  schedule. 

Consider  a  certain  component  of  an  operating  system. 

The  component  can  fail  at  a  random  time  Y;  occasional  inspection 
is  necessary  to  determine  when  it  does  fail.  The  cost  of  a 
failure,  vt,  where  t  is  the  time  between  the  failure  and  its 
inspection,  and  v  is  the  cost  per  unit  of  t,  can  be  reduced 
by  frequent  inspection.  However,  there  is  a  cost  c  for  each 
inspection,  so  the  number  of  inspections  must  be  kept  small. 

The  best  compromise  between  these  two  conflicting  requirement 
in  the  sense  of  minimizing  expected  total  costs  is  called  the 
optimum  inspection  schedule.  In  a  finite  horizon  problem  we 


assume  the  system  will  stop  operating  when  a  stipulated  time 
T  is  reached  or  when  a  failed  component  is  detected. 

In  this  section,  we  will  derive  a  formula  for  supremum 
expected  total  cost  for  any  inspection  schedule  assuming  the 
mean  of  Y  is  given. 

Definitions: 

U  is  the  mean  of  Y; 

is  the  class  of  all  distribution  functions 

I 

for  nonnegative  random  variables. 

O',,  is  the  subset  of  that  have  expected  lifetime  n; 

•  H*  I 


p  2  r 

; »  IS  the  subset  of 


increase ; 


that  have  exactly  two  points  of 


X  =  (Xq,x^,.,.,  x^,  denotes  an  inspection  schedule 

(short:  schedule)  ,  where  Xq=0,  =  T,  and  x^  < 

l”0,l,2, . . . ,n7 

m  is  the  index  of  an  inspection  in  x  such  that  x^  <  y , 
’^m+l  ^  ^  ' 

C(x,F)  is  the  expected  total  cost  of  x  if  F  is  the  c.d.f. 

of  Y; 


The  supremum  expected  cost  of  x  is  u(x)  =  sup  C{x,F). 

Fet| 

If  <^6J^otes  two  increasing  points  of  the  c.d.f., 

2 

F(yi,y2),  of  Y,  which  is  a  member  of  then  the  probability 
function  of  Y  is 


P(Y  = 

P(Y  =  y2) 


^2 


-  y 


y 


1 


-5 


9 


< 


where  0  <  <  p  < 

Proposition  2.1.  If  x  is  a  given  schedule  with  x  <  ' 

m+l  — 

and  F(y^,y2)  is  a  c.d.f.  in  ^  with  V  <  )l 2  -  ^ni+1 
C  (x,F  (y^^  jy^)  )  <  u  (x)  . 

Proof. 


C(x,F(y^,y2)  ^ 

=  C  t  V(X.^J  -  +  Km+DC  yj)) 

,  if  Xi  <  yi  <  and  u  <  y2  1 

By  taking  partial  derivative  with  respect  to  y^,  we  have 

C(x,F(y,  ,y  )  ) 

3y2  ^  ^ 


which  shows  that  C(x,F(y^,y2))  ®  strictly  decreasing 

function  on  y^;  hence,  we  have 

C(x,F(y^,y2)) 

<  C(x,F(y^,U+)) 

=  (ra+1)  C  +  V  -  ^^) 

<  (m+l)  C  +  V  (x^^,  -  X  ) 

m+l  m 

=  lim  C (x,F(x  +,  y  ) ) 

\/  -von  ^ 


<  U(x) 


Q  •  G  •  D  • 


Proposition  2.2.  U(x)  ^  sup  C(x,F) 

^45 


Proof.  Let  (!  *  =  F  is  a  step  function}  and  d(F,G)  = 

suplF(t)  -  G(t)l;  theorem  2.2.  in  Hoeffding  [3]  states 
*t 

(2.1.1.)  sup  C(x,F)  =  sup  C(x,F)  =  sup  C(X,F) 

^2  y*  Fe(^ 

if  (A)  C(x,F)  is  a  continous  function  on(^^  for  x  fixed,  in  the 
sense  of  d(.,.),  and  (B)  for  any  Fe^^,  there  is  a  sequence  {F^}^C^ 

such  that  lim  =  F  in  the  sense  of  d(.,.). 

n-^oo 

The  satisfaction  of  A  and  B  are  easily  verified. 

Definitions : 

=  C(x,F(y^,y2) ) 

y-  -  P  V  -  Yi 

=  %‘yi'  *  9j<y2' 

,  where  x  <  y  <  x_. . ,  1  <  i  <  m,  m+1  £  j  £  n+1, 

'  m  —  m+i  “  ~ 

gj.(y)  =  (r+l)C  +  V(x^^3^  -  y)  if  y  e  (x,  and 

r  ^  0,1, ...,n 


=  (n+l)C 


if  y  >  T  and  r  =  n+1  ; 


Gi.(x)  =  f^.Cx.  +,  X.+)  ; 

Gj.(x)  =  * 

Theorem  2.1.  For  any  schedule  x, 

u(x)  =  sup  (G^ j  (x)  ,Gj^  (x) } . 
0<i£m 

m+l< j<n 
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Proof.  The  partial  derivatives  with  respect  to  and 
on  respectively  are 


°^n‘Yi-y2> 


y,-^ 

(2.1.3)  =,-i —  -[(i-j)C  +  V(x...-x  )]<  0,  if  j  <  n, 

iy2  1+1  3+1 


y2-M 

(2.1.4)  [(i-n)  C  +  V(x^^j^  -  y2)]<  0,  if  j  =  n+1 . 


^^i<yi>y2- 

3y2 


u-y, 

(2.1.5)  ={y~ry^)2  [(i-j)C  +  +  0,  if  j  ^  n, 

(2.1.6)  ~  (y2~yx^  ^  [(i-n)C  +  V(x^^j^  -  y^^)  ] ,  if  j  =  n+1. 


If  the  value  of  (2.1.6)  is  negative  then 
^in+l^^l'^Z^  <  '  together  with  (2.1.1),  proposition  2.1., 

(2.1.3),  (2.1.4),  (2.1.5) ,  we  have 


U(x)  =  sup  (G. . (x) } 

O^i^m 

m+l^j^n 

If  the  value  of  (2.1.6)  is  nonnegative  then,  we  have 

n+l'Yl'Y2’  -  'i  n+l'Yl-y2'  =  ’i'Yl*  i  ^i'*’ 


which  implies 

U(x)  =  sup  {G. .(x),  G. (x) } 
0<i<m  ^ 

m+l£j^n 
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Q.E.D. 


iMmM 


Utioiiiiafe 


2.2.  The  minimax  schedule  is  a  schedule  with  equally  average 
increments 


Definitions : 

i  denotes  the  index  of  an  inspection  of  a  schedule  x  such 
that  0  ^  i  ^  m; 

j  denotes  the  index  of  an  inspection  of  a  schedule  x  such 
that  m+1  <  j  ^  n . 


Proposition  2. 
j=m+l  but  ij^m, 
G,  . (x)  is 
function  on  x. 

X 


3.  X  is  a  given  schedule.  If  j  7<^  m+1,  or 
then 

an  increasing  function  on  +  ^nd  decreasing 


If  j  =  m+1  and  i=m,  then 

^m(m+l;  ^  decreasing  function  on  x. 

Proof . 


If  j  ^  m+1,  or 


9G.  .  X  .-VI 
=  -2 


3x . 
1 


(X .-X. ) 

3  1 


j  =  m+1  but  iji  m  then 
{(i-j)C  -  V(x.^3^-x^^j^}  <  0 


3G^ 


3x 


i+1 


0 


3G. 


U  -  X. 


iJ_  =  V  - - -  -  >  0. 


^^j+1 


X  .  -X  . 

3  1 


If  j  =  m+1  and  i  =  m  then 
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8G  (m+l) 
m 


m 


,2  f-C  -  V(x^^2  -  <  0 


(X 


Itt+1 


m 


Q.E.D. 


Proposition  2.4.  For  any  schedule  x,  ^ (x)  <  G^ (x)  if  and 
only  if  Gj  (x)  <  Gj^(x). 

G.  (x)-G. (X) 

Proof.  G.  .(x)  =  G.  (x)  +  - - -  (u  -  X.) 

<  G^(x) 

Gj  (x)  <  G^(x),  since  y-x^  >  0 

Q.E.D. 


Proposition  2.5.  If  x  is  the  minimax  schedule  and  (j)  n  >  j  ^  m+l 

is  nonempty  set,  then  u(x)  =  sup  {G. . (x) } . 

i,j 


Proof.  If  u(x)  =  G.  (x)  for  some  in<Tn,  and  G.  (x) 

^0  0 

>  G.  . (x)  for  all  j  e  {j}  n  >  j  ^  m+l,  then,  by  nonemptiness 

of  {j}  n  ^  j  ^  m+l  and  proposition  2.4.,  we  have  G  ,  (x)  <  G.  (x) . 

m+i 

Hence,  we  can  find  a  sufficiently  small  e  >  0  and  the  corresponding 
schedule 


x^  =  (0  , . . . ,  x.-  ^ 


i"2  (m-i+1)  '  •  •  • '  ’^m+l”^'  ’‘m+2,...,  x  )  such 


n 


that 

(2.2.2)  G^^^(x)  <  <  G.^(x). 

Since  ^i+1  “  ^i'  i=0,l,...,  m  and  ■■  ^  ^j+l'^'j' 
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j  =  m+2,...,n,  we  also  have 

(2.2.3)  G.  (x*^)  <  G.  (x)  £  G.  -  (x)  ,  i=0,l,. . .  ,in  and 

(2.2.4)  G. (x^)  =  G. (x)  <  G.  (x) ,  j  =  m+2,...,n. 

J  ^  0 

Case  1.  X  . ,  >  p . 

m+i 

If  e  is  sufficiently  small  such  that  -  e  >  y,  then, 

by  (2.2.2),  (2.2.3),  (2.2.4),  proposition  2.4.  and  theorem 
2.1,  we  have 

(2.2.5)  U(x‘-)  =  suplG..  (x^),  G^(x^)} 

i/ j  ^ 

<  G.  (x)  =  U  (x)  . 

^0 

Case  2.  =  V* 

In  this  case,  -  £  <  y.  (2.2.5.)  becomes 

U(x^)  =  sup  {G.  .(x^),  G.  (x'^)} 

O^i^m+1  ^ 

m+2^j^n 

<  G.  (x) 

0 

=  U (x)  . 


Both  case  1  and  case  2  lead  to  a  contradiction;  thus, 

U (x)  =  sup  {G. . (x) } . 
if  j  ^ 


Q.E.D. 


Definitions : 

If  X  is  a  schedule,  G^ (x)  - 
with  respect  to  (i,j)  and  Ax(j,i) 


G^ (x)  is  called  an  increment 
G. (X)  -  G.(x) 

_  - i -  ig  called 

Xj  -  x^ 
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an  average  increment  with  respect  to  (i,j).  If  x  is  a 

schedule  such  that  for  some  Z  Ax(r+l,r)  =  Z,  r=0,l,...,n 

then  X  is  called  a  schedule  with  equal  average  increment  (SEAI) . 

Proposition  2.6.  If  for  a  schedule  x,  there  exists  i^  and  jp 
such  that 

(2.2.6)  G.  1  (x)  =  ...  =  t?.  „(x)  and 

iQin+l  ign 

, 2.2.7,  .  ... 

then  X  is  a  SEAI . 


Proof.  For  convenience,  let  =  G^^ (x) ,  G^  =  G^(x), 

A(i,j)  =  Ax(i,j).  Since 


^ij  “  (U-x^) 

=  Gj  -  A(j,i)  (Xj-y) 

for  any  schedule  x  and  any  existing  pair  (i,j)  we  have, 
by  (2.2.6;  , 

A(j,ig)  -  A(j+1,  ip) 


G. 


-  '''^0 
u  -  X, 


y  -  X. 

A 


=  0 


,  3  “  m^"  1,...,  n. 

If  jjj  f  m+1  or  Xj^^^  >  y,  then,  by  (2.2.7)  and  the  same 
reason  as  above,  we  have 

(2.2.8)  A(jQ,i)  -  A(jQ,  i+1)  =  0  ,  i=0 , 1 , . . . ,m-l .  If  jQ=m+l 
and  *  y  then  i*0,...,m,  which  implies 

*  ®ni+l  '  (2.2.8)  still  holds. 
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In  other  words,  the  following  holds; 

(2.2.9)  A(m+1,  Iq)  =  ...=  A(n,iQ)  and 

(2.2.10)  A(jQ,0)  =  ...  ==  A(jQ,m). 

Since  (2.2.9)  and  (2.2.10)  have  the  same  terra  A(jQ,iQ), 
we  have 

(2.2.11)  A(jQ,0)  =  A(jg,l)  =  ...  =  A(m+l,iQ)  =...  =  A(n,iQ). 

Setting  the  expression  of  (2.2.11)  equal  to  Z,  by  the 
definition  of  A(.,.)  and  A(jQ,i)  =  A(jQ,i+l)=Z,  we  have 


G..  “  G.  —  G.  “  G.  +  G.,1  “  G. 
1+1  1  3.  1  1+1  3. 


=  (x.^  -  x,)Z  .  (x,^,  -  x^^)  2 


=  (x^^l  -  X.)2. 


This  implies 

(2.2. 12. a)  A(i+l,i)  =  Z  for  i=0,l, . . . ,m-l;  similarly,  by 
A(j+1,  ig)  =  A(j,ig)  =  Z  we  have 

(2.2.12.b)  A(j+l,j)  =  Z  for  j=m+l,...,n. 

What  remains  to  be  shown  is  that  the  value  of  A(m+l,m) 
is  also  Z.  This  follows  from 


Z  = 


G.  -  G. 


^0  -  ’‘i. 


=  (G.  -G  .i)+(G.  -  G,  • 


.  . +  (G  . 1 “G  ) 
m+1  m 


-*i 


+(G  -G„  .) 
m  m— 1 


“  G,  ) 
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Z  (x  •  “  X  .  1  +  X  .  1  ~  X  .  -  »  .  .  ^—4.1  )  ^ 

jQ  ^0  ^  iTi+l  m+l  m 


X  .  -  X. 

3n 


+  Z  (x  -x_  , +  ...  -  X . 
m  m-i  1, 


Z(x,  -  x^ . , )  +  (G  ,  -  G  )  +  Z{x  -X. ) 

m+l  m+l  m  m  i_ 


X  .  -  X  . 

^0  ^0 


=  z  + 


G  .  ,  -  G  -  Z  (x  -  X  ) 
m+l _ m _ m+l _ m_ 

X  .  -  X  . 

In 


'0  "0 

Cancelling  Z  from  both  sides  of  above  equation  we  have 

G  . ,  -  G  =  Z  (x  ,  ,  -  X  ) 
m+l  m  m+l  m 


which  implies  A(m+l,m)  =  Z.  Together  with  (2. 2. 12. a)  and 
(2.2.12.b)  the  proposition  is  proved. 


Lemma  2.1. 

Suppose 

a,b,A,B  are  real 

numbers  and  B  +  b  >  0 

B  >  0 ;  we 

have 

a)  b  >  0; 

A  + 

< 

A  if  and  only  if 

a  -  A 

B  + 

b 

B 

b  B* 

b)  b  <  0; 

A  + 
B  + 

■!< 

A  if  and  only  if 
B 

a  .  A. 

F  ^  B 

1 

\ 

\ 

j 
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Proof. 


rH  ‘  I  ‘ 


<=>  a  <  ^ 
B 


«=>  I  <  I  (by  b  >  0) 


^  ^  <  b.  ^ 

B  +  b  B^~'  B 


<=>|  >  g  (by  b  <  0) 


Definitions : 


'4) .  (x)  =  sup  {G,  .  (x)  } 

J  f\  ^  i  J 


0<i<m 


Q.E.D. 


$ . (x)  =  sup  {G. , (x) } . 

^  m+l^j<n 

Theorem  2.2.  A  minimax  schedule  x  is  a  SEAT  with  some  non¬ 
negative  increment  Z. 

Proof.  We  prove  the  theorem  under  the  following  two  cases; 

Case  1;  If  ^j^n>j>m+l  ~  ^  then  x  is  a  SEAl  with  average 
increment  0. 

Case  2;  If  /  <(>  then  x  is  a  SEAI  with  some 

n>3^m+l 

nonnegative  average  increment  Z. 

Case  1.  By  theorem  2.1.  and  =  <!'»  we  have 

U(x)  =  sup  {g. (X) }. 

0<^i^n  ^ 

Suppose  i  is  the  largest  index  such  that  U(x)  =  G.  (x) 

* 
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and  suppose  i  ^  n.  Because  G. (x)  is  a  continuous  function 

X  X 

in  *i+l'  increasing  in  decreasing  in  x^  together 

with  the  assumption  above,  we  can  find  a  sufficiently  small 

e  >  0  and  a  schedule 

X  =  (x.fX-,***,  X-  ,x.  -  e,...,x  ) 

^  ^  ^x  x+1  " 


such  that 


G.  (x^^)  <  G.  (x)  =  U(x)  still  holds,  and, 
^x+1  ^x 


(2.2.13) 

G. (x^)  <  G  (x)  =  U(x) . 

^x  ^x 

Thus  U(x^)  £  U(x)  ;  however  x  is  minimax  so  U(x^)  =  U(x) 
and  ij^e  <  i^^. 

Repeating  the  above  procedure  at  most  i^^  +  1  times, 
we  reach  a  schedule  x*^  such  that  U(x*^)  =  U(x)  and  i^®  =  0. 

Again,  let  n  0  be  small  enough  with  the  corresponding  schedule 

X  =  (0,  Xj^-n,  X2,...,x^) 

*  0  *  0 
such  that  Gj^(x  )  <  Gq  (x  )  and  Gg  (x  )  <  Gg  (x  )  then  we  have 

U(x*)  <  U(x*^)  =  U(x),  a  contradiction.  Thus, 

If  there  is  a  i-  such  that  G.  (x)  <  G  (x)  =  U(x) 

0  ig  n 

then  we  can  find  a  sufficiently  small  e  >  0  along  with 
a  schedule 


such  that  G. 

^0 


and 


X  -  (x,,X2,...»  X.  ,x.  .  +  e,x.  ^2 

1  ^  ig  ig+i  ig-^-: 


(x  )  <  u(x)  still  holds.  Since 


Gj,(x  )  =  Gj,(x),  r=0,l, . . .  ,ig-l,ig+2, . 


.  ,x  ) 
'  n 


n. 
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X  is  also  a  minimax  schedule.  By  induction,  we  can  find  a 

*  *  * 

minimax  schedule  x  such  that  G  (x  )  <  u(x  )  which  is  a  con 

n 

tradiction.  Hence,  G^(x)  =  u(x)  for  all  i. 


Case  2. 


Case  2.  can  be  proved  by  three  steps: 

Step  1.  To  show  ({),.(x)  =  ...=  ({)  (x) 

u  m 

and  ...  =  4»^(x). 

Step  2.  To  show  there  exists  i^  and  such 

that  G.  ,  (x)  =  ...  =  G.  (x)  and  G,.  (x)=. 

loHitl  iljn  0.^ 

Step  3.  To  show  Z  is  nonnegative. 


Step  1.  By  proposition  2.5.  0(x)  =  sup  {G. . (x) } 

i»'j 

and  U(x)  =  sup  {i(».  (x)}  follows.  By  proposition  2.3. 
in+l^j<^n  J 

and  the  definition  of  4/^  (x)  ,  we  know  ^  (x)  is  continuous 
and  decreasing  on  x^  for  j  =  m+2,...,n  and  is  continuous 
and  increasing  on  x^^^^  for  j=m+l,...,  n-1.  If  we  treat 

as  a  fixed  initial  point  then  the  method  in  case  1  can 


show 


"^m+l^^^  =  ...  =  =  u(x)  . 


;  analogously,  if  we  treat  x  ,  as  a  fixed  terminal  point 


we  have 


<t)Q(x)  =  ...  =  4>j^(x)  =  u(x). 


step  2.  Suppose  ^^j(x)  and  G^  ^^j(x) 


Since 


=i.m+2<*>  =  ■*  -  ’‘i  >  ^  X  X. -  “‘■’‘in’’ 

0  0  0  in+2  In  0 


G. 


i-m+Kx)  =  (i  +1)C  +  V{x.  .  -  x.  )  +  - - — 

0  O  +  ^  .  1  “  X  . 

0  0  m+1  1 


G  .  ,  (x)-G. 

^  ^0  (y-x.  ) 

^0 


0 


we  have 


V2'*'  -  °i  W"’  -  ®1  <*> 
<■>  - 2 -  <  0 


X  . «  -  X. 

m+2  1, 


X  ^ ,  -  X. 

m+l  1, 


Since  x^^2  ~  ^m+l  ^  ®  lemma  2(a)  ,  the  above  inequality 
is  equivalent  to 


(2.2.U)  G^„(x)  -  C.^^j(x)  -  G  lx) 

- -  -  I  <  - - 


%+2  -  ^m+1 


X  , ,  -  X  .  . 

m+l  ip 


At  the  same  time,  by  the  definition  of  have 

G.  , (x)  <  G.  (x)  for  all  i,  and  since 
im+1  —  IqHi+i 

=lm*l<=‘'  =  -  X.  " -  ■ 

m+l  1 
G  ,  (x)  -  G.  (x) 

G..  -xi  (x)  =  G_.,{x)  ^  ^0 

°  *m+l  "'ip 

we  have 


(2.2.15) 


^m+1  “  ’'i  ^ro+1  ■  ^i„ 
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Equations  (2.2.16),  (2.2.17),  and  proposition  2.6.  implies 
X  is  a  SEAI. 

Step.  3.  If2<0,  byxisa  SEAI,  we  then  have 
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G. (x)  -  G. (x) 

G..(x)  =  G.(x)  +  ^  _  ;r. - (v-i) 


=  G^ (x)  +  Z (u-x^) 


<  G^ (x)  for  all  (i,  j)  , 

which  contradicts  u(x)  =  sup{G. .(x)}. 

Thus,  the  theorem  has  been  proved. 

Theorem  2.3,  If  x  is  a  SEAI  with  non-negative  average  increment 
Z  then  Gq  (x)  +  Z|i=u{x)  . 

Proof:  If  Z=0,  by  the  definition  of  SEAI,  we  have 

Gq  (x)  =Gj^  (x)  =  ...»G^(x}  which  implies 

(x)=G^(x)=Gj (x)=Gq (x)  for  all  (i,j);  thus  by 

theorem  2.1., 

U(x)=  sup{G^j(x),  G^{x)}  =  Gq(x) 

^  f  j 

=  Gp(x)  +  Zu. 

If  Z  >  0,  by  definition  of  SEAI,  we  have 
Gq(x)  <  Gj^(x)  <  ...  <  G^(x)  which  implies  ^  (x)  >  G^{x) 

for  all  (i,j);  thus,  by  theorem  2.1,  we  have 

u(x)  =  sup{G..(x),  G. (x)}  =  sup{G..(x)}. 
i»j  ^  ifj  ^ 

Because 

Gij(x)  =  G^{x)  +  A(j,i)(M-x^) 

=  Gj (x)  -  &(j,i) (Xj-p)  for  all  (i^j); 
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G. (X)  -  G^{X) 

(2.2.18)  A(j,i)  =  -J — tie - 

Xj  X. 


■  (x)-G._^  (x)+G^_3^(x)-...+G.^^^(x)-Gj^(x) 


^j-1  '^j-2  +•••+ 


^  ‘i-2  ■^••■■^^1  ^ 

5.  ,  +6.  _+...+6. 

D-1  D-2  1 


the  following  hold 

(2.2.19)  G.j^^j^(x)  =  =  ...  -  G.^(x)  for  all  i 

(2.2.20)  GQj(x)  =  G^j(x)  =  ...  =  Gj^^(x)  for  all  j. 

Equation  (2.2.19)  and  (2.2.20)  imply  all  G^^j{x)'s  are 
equal.  This  is  evident  by  the  following  diagram. 


”0  r.+l 


S  b+2 


••  *  Sn’** 


■^1  ni+l'’'’ 


•’l  n.+2  '*> 


'=m  '  '’n,  ™.2<''>  '  •••  - 

,  where  the  row  equalities  follow  by  (2.2.19)  and  the 
column  equalities  follow  by  (2.2.20) . 

By  (2.2.18) ,  we  have 

U(x)=sup{G.  .  (X)  }=  G  (X)  =  Gq(x)  +  A(n,0)u 


=  Gq(x)  +  Zy 

Q  •  E  •  D  • 
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2.3.  The  properties  of  objective  functions. 


Let 

g(b,n) 


c+ 


bTv  -  (n+Dc 
(l+b)’’'^^  -  1 


+  ^  +  bvp 


*  * 

Proposition  2.7.  Ifb  >0,n  is  a  positive  integer,  and 

*  * 

g(b  ,n  )  =  min  g(b,n),  where 
(b,n)  eA 

A  ={(b,n)  IbcR^,  nfr|l,  and 

(bn  -  b^  -  1)  (1+b)"  +  1  <  0} 


then 


★  ★ 

X  -  b  T  -  (n  +  l)v  c 

°0  *  n*+l  *  ' 

^  (1+b  -  1  b  V 


=  (l+b*)""  6  -  ^  £  , 

u  b 


* 

for  r=0,l,...,n  ,  and  define  the  minimax 

schedule. 

★ 

Proof.  By  theorem  2.2.  and  2.3.,  the  minimax  schedule,  x  , 
minimizes 

Gg  (x)  +  Zxi 

where 

-  G^(x)  =  Z6^ 

6^  >  0,  Z  >  0,  r  =  0,l,...,n. 


and 


"r-H 


n 

I  ^=T 
r=0  ^ 


Equivalently, 

(2.3.1)  C  +  viSg  +  Zu 

Is  minimized  over  6g  with  constraints 
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(2.3.2)  C  +  =  vfij.  +  Z5j.  =  (v+Z)  6^. 


(2.3.3)  n 

I  6  >  0 ,  Z  >  0 ,  r=0 , 1 , . . . ,  n 

r=0  ^ 


holding. 

However  (2.3.2)  and  (2.3.3) 
(2.3.4) 

and 


n 


Z5=T  6  >0,Z>0,  r=0,l,...,n. 

r=0  ^  ^ 


Let  b  =  then  from  (2.3.4)  we  have 

V 

(2.3.5)  5^=(l+b)^  (Sq  £  ,  r=0,l,. 


n 

E  6=T,6  >0,b>0 

r=0  ^  ^ 


which  imolies 
n 


-  ^  n  iv,  1  ^  ^  (1+b)  ^  c 

5  E  (1+b)  -  1—5 - 

'"r=0  r=0 


-  =  T, 


V 


6  >  0,  b  >  0. 

r 


Rearranging, 

‘«0  -  !?>  T.  '  «r 

r=0 

and  solving  for  6^^  we  have 


(2.3.6)  «o  •  5?  *  >■  ^ 


n 

E  (l+b) 
r=0 


=  £_  +  (T  -  - 

bv  “TT— 
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bT  -  (n+l)£  ^ 

-  1  * 


6  >  0 ,  b  >  0 . 

r  ' 


From  (2.3.4)  in  order  for  ^  0 ,  we  have 

(1  +  -)  6  -  £  >  0;  i.e.  6^  >  §/ (1+|)  >  0. 

V  r  V  r  V  V 

Thus,  if  6  >0  then  6.,  5,,...,  6  >  0.  From  (2.3.5) 

n  u  i  n 

and  (2.3.6)  we  get 

6n  =  (1+b)"  fig  -  ~  -  I  ;  hence, 


6  >  0 
n 


<=^ 


(1+b)"  6q  > 


(l+b)*^  -1 
b 


c 

V 


<=^ 


>  -  1 
b(l+b)" 


c 

V 


£  n 

bT  -  (n+ljv  ^  c  ^  (1+b)  -  1 

(l*b)"*l  -  1  bu+b)" 


C 

V 


<=>  (nb  -  b^  -  1)  (1+b)"  +  1  <  0. 

Thus,  we  get 

c  +  =  (v  +  Z)  5^ 


n 


I  5  =  T, 
r=0  ^ 

(i.e.  (2.3.2) 

6> 


5 


0 


bT  ~ 
(i+b) 


5^  >  0,  Z  >  0,  r=0,l 


and  (2.3.3)  holds) 


c 

(ntl) V  c 

n+l_  ^  bv  ' 


n 


(nb  -  b^-1)  (1+b)"  +  1  <  0, 
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g(b,n')  -  g(n'-l)  >  0 

from  which  the  proposition  follows. 

Since  (b,n)  e  A  if  and  only  if 

(nb  ~  b^  -  1)  (1+b)”  +  1  <  0,  we  define 

c 

fn(b)  =  (nb  -  ^  b^  -  1)  (l+b)*^  +  1.  We  then  have 

n(b)  =  max{n|f  (b)  <  0}. 
neljl  " 


Let  n.  =  lim  n(b). 
^  bH-0 


0  r 

Proposition  2.9.  n.  >  max  {n|n(n+l)  <  — . 


Proof.  Expanding  (1+b)  we  have 

f  (b)  =  .  11  +  o(b^)  <  0. 

n  2  c 


Then 


n(n+l)  Tv  .  0(b  )  . 

— 5 —  -  —  <  —  ■  -  ,  as  b->^0  and  thus 


n(n+l)  £ 


Q.E.D. 


The  equality  in  Proposition  2.9.  holds  when  2Tv  is  not 

2Tv  ^ 

an  integer,  and  max  {n|n(n+l)<  — is  the  optimal  number  of 

neljl  ^ 


inspections  in  Derman's  schedule,  the  minimax  inspection 
schedule  when  nothing  is  assumed  about  the  distribution  of  Y. 

2Tv 

Proposition  2.10.  lim  f  (b)  =  n(n+l)  -  — =— . 

.  n  c 

b  -0 

Proof.  Follows  from  the  proof  of  proposition  2.9. 

Q.E.D. 

From  the  definition  of  Hq  and  proposition  2.10.,  we  have 
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as  follows 


f  (0  )  >  0;  thus  we  can  define  b 

”o+i  "o+i 

Definition:  b  =  sup{blf  (b)  <  0,  f  (b)  >  0}. 

■^0+1  b>0  ''o  0+1 

Proposition  2.11.  b  is  the  unique  zero  of  +  (•) 

"o+l  "o-rl 

on  (0,«>)  ,  and  n(b)  =  n-  for  0<b<b 

^  ~  0+1 

Proof . 

(2.3.8)  f^^  (b)  =  b(l+b)"0[  (n^  +  l)  (np+2)  -  {n^+3)'^] . 

t 

By  (2.3.8)  we  have  f  +  b  and  f  (oo)=-c«  which  implies 

no+1  n^+l 

f  (•)  has  an  unique  zero  on  (0,‘»)(see  Fig  1). 

"o+i 

Let  f  (x)  =  0,  we 
^0+1 


have 


Fig.l. 


(2.3.9)  '  (b)  >  0  if  0  <  b  <  X  r ^ 

“^O+l 

(2.3.10)  f  (b)  <  0  if  b  >  X. 

0+1 

Also,  by  proposition  2.9.,  we  have 

(2.3.11)  =  b(l+b)"0[nQ(nQ^j^)  -  ^  ° 

which  implies 

(2.3.12)  f  (b)  <  0,  if  0<b<x. 

"o 

The  definition  of  n(b),  b  ,  (2 . 3.9 ) , (2 . 3 . 10) , 

■^O+i 

(2.3.12)  prove  the  proposition. 

“  =  "•  ‘k+l**’  <  0  for  X  >  and 

fj^+l(x)  >  0  for  0<x<bj^^j^,  then  we  have 

*kf2<'’kfl  *  •  '("^I'bk+l  ♦  ♦  ^k+l  ♦  ^  ">kn 


+c)^  -  1}  x  (1+b^+l  +  +  1-  AS  e-0 
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(2.3.13)  t^^2(b^^l+)  i  0 

holds;  also, 

We  can  inductively  define  k=nQ^j^,...,  as  follows: 

Definition;  bj^+2  ~  I  (b)  <0  »  ")c+2  ®  ^  * 


Proposition  2.12.  The  unique  zeros  b  , . . .  of  f  ,..., 

■^0+1  0+1 

indicate 

n(b)=k  if  bj^  <  b  £  where  k=nQ  ,nQ^j^ , . .  and  b^  =0. 

Proof.  f^+2(b)  =  b(l+b)’''^^[{k+2)  (k+3)-^  -  (k+4)2!^], 

(2.3.13)  imolies  bas  a  unique  zero 

on  (bj^^j^,«>)  .  Denote  the  zero  by  x  ;  then  we  have 

fk+2^‘^^-‘^  if  bj^^j^<b<x, 

^k+2  b>x. 

Also,  from 

f^^l(b)  <  0  for  b  >  b^^^ 
together  with 


^k+1 ^^k+1^ 


0 


we  have 

f^^j(b)  <  0,  <  b  <  X. 


Q.E.D. 


Proposition  2.13.  g{b,n(b))  is  a  continous  function 


Proof 


lim 


g (b,n(b) ) 
k+1 


lim 

b-*-b 


g(b,n(b)  ) 
k+l 


=  *'k-n"’k+l' 


irn" 


=  0 


-  1] 


[(1+b,^)'"  -  1] 


Q  •  E  •  D  • 


Proposition  2.14.  For  a  given  b>0,  let  k-n(b)  and 
k+1 

s(b)  =  —  /  if  Tv  s(b.  , )  >  6  then  g(b,k) 

(l+b)*^  ^  c 


is  a  convex  function  on  1  * 

Proof.  By  taking  first  and  second  derivatives  on  g(b,k) 

with  resoect  to  b,  we  have 

(2.3.14)  d  g(b,k) 

3F 


'Tjd  +(l+b)’^  {kb  -  ^ 


(k+l  2 


-  1} 


=  ( 


[ (1+b) 


k+1 


-  11' 


c 


and 

(2.3.15) 


— ,  g(b,k) 
db'^ 


~Tl  (i+b)^"^ 
=  {_S - 


{ (l+b)^"^^  4)  (b,k)-  (k^+3k+l)  (l+b)-k(k+l)  (1+ 
{(l+b)^"^^  -  1}^ 


c 


f  where  <p{b,k)  =  k(k-l)b+2 


(k+1)^ 

(k-2) 

c 
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Let  X  =  by  definition,  we  have 

1+  (l+x)^'^^{  (k+l)x  -  x^  -  1}=  0 


^  .  .  S|0 


Tv 

For  convenience,  let  a  =—  and  s  =  s{x),  then 

-x^«^(x,k)  =  (sa-2)ax^  -  ( 3s+2)  ax^+2as^x^-3s^x 

+s^  =  a{(sa-2)x^  -  ( 3s+2)  x+-^s^ }x^+s^ ax^-3x+s}. 

The  discriminant  of  the  Quadratic  polynomial  in  the  first 
brace  of  the  above  formulation  is 

(2.3«16)  (3s+2)  —  4(sa—2)'j-g-  <  25  +  — *■  sa 

,  and  in  the  second  brace  is 
(2.3.17)  9-||as. 

The  condition  as  ^  6  makes  both  (2.3.16)  and  (2.3.17) 
negative;  hence,  we  have  <|)(x,k)  <  0.  Since  <(>(•, k)  +  b 
for  <  b  <  =  X,  we  then  have  £  g(b,t)>  0  for 


db 


Q  •  E  •  D  • 


Proposition  2.15.  g,b,„(b))  >0  if  b  >  /? 


Proof,  b  >  /g~ 
/  uv 


(2.3.18)  -  iy  +  XE  >  0 

b  ^ 


Also, 

c 
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-  1  +  n+  f  (n-H)  Tv  2  _  ,, 

-  J-  '  Tv  ^  Tv  ”  c  '  Tv^  ” 


1  +  <  0 


which  implies  n+1  ^  ^ 

■  wi  >  b  .  %  >  a 

Tv  n+1 


;  thus,  l+(l+b)"{nb  - 


-  1} 


<  l+(l+b)”{nb  -  -  1} 

c 

=  fjj(b)  <  0,  where  n=n(b). 

Together  with  (2.3.18)  and  (2.3.14),  we  have 
^  g(b,n(b)  )  >  0. 

O.E.D. 


Proposition  2.16.  Tv  „  ,k\  a 

r—  s  (b)  t,  . 
c  b 


Proof.  Recall  that  s(b)  = 


H+T"  ' 


n  =  n(b).  By  proposition  2.12.,  we  have  n(b) +j^  and 


^n+1  ''^b* 


Q.E.D. 


As  a  summary  of  proposition  2.13,  2.14,  2.15  and  2.16, 
we  have  the  following  theorem. 

Tv 

Theorem  2.4.  If  s(0+)  >  6,  then  the  objective  function 
g(b,n(b))  is  continous  and  piecewise  convex,  and  the  optimal 

*  -..I.. 

b  is  in  [  0 ,  /  c  ) . 

/  w 
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Proof.  By  proposition  2 . 13 . ,2 . 14 . , 2 . 15 . ,  2.16.. 

Q.E.D. 

2.4.  The  asymtotic  minimax  schedule  as 
Definitions:  denotes  a  schedule  on  [0,T]. 

For  a  given  schedule  x,  n(x)  is  the  total  number  of  inspections 

assigned,  6^ (x)  is  the  length  between  ith  and  i+1  th  inspections, 

and  u(x)  is  the  corresnonding  supremum  expected  cost.  If  x 

is  a  SEAI  then  Z(x)  is  its  average  increment  level.  If 

6.  (x)=  ,  I  ,  ,  for  all  i,  then  x  is  called  a  strictly  periodic 
1  nTx)  +1 

schedule.  If  s  is  a  real  number,  [s]  denotes  the  largest 
integer  which  is  small  then  s. 


Proposition  2.17. 
then  X  is  a  SEAI. 

Proof.  A(i+l,i)= 


If  X  is  a  strictly  periodic  schedule 

(i+2)c  +  v6^^^ (x) - (i+1) c  -  v6^(x) 

(x) 


c(n(x)+l) 

T 


Q.E.D. 


Proposition  2.18.  If  {T^ }  is  a  time  sequence  and 

*  * 

as  then  n(x,p  ) -♦«’  as  where  is  the  minxmax 


schedule  on  [0,Tj]. 


Proof.  Suppose  x^  is  Derman's  schedule  on  [0,?^] 


then  Z  (x)p  )  =  0,  by  proposition  2.12.,  we  have 

0  * 

n(x„  )  <  n(x.j,  )  . 

j  1 
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But  lim  n(x  ^  )  =  «>.  The  oroDosition  follows. 

3 

Theorem  2.5.  If  is  the  minimax  schedule  on  [0 
there  exists  a  strictly  periodic  schedule  x^  for 


that 


lim  u(x^,)  =  lim  u(x!,) 

T-*-®'-'  T-*«>  ^ 


-  C  +  V.'MC  u  /JJc  . 

/  ~  ^  /“ 


Proof.  By  theorem  2.3.,  x^  is  a  SEAI  and 
(2.4.1)  u(x*)  =  c  +  v^q(x*)  +  Z(x*)u. 

Let  {T  }  be  any  time  sequence  such  that  T. -*•<*' 
^  3 

and  x°  be  a  strictly  neriodic  schedule  such  that 
j 


1  +  n(x°  )  =  j 
3 


/JJC 
/  V  j 


u(x^  )  =  c  +  V 


T . 


/P 


j  Tj  ! 


1  i 


+  U 


^  J 


T 


Since 


r  ^ 


1i« 

im 


|/  V  J 


0 


=  .'-r^  and  u  (x  )  1  u  (x  )  , 
3-*"  T.  /!!£ 


,T1  then 
each  T  such 


as  3 


we  obtain 


(2.4.2)  lim  u(x„  )  <  c  +  Vy  |ic  +  p  ^ 


T . '  - 
1 


7ir 


/I 


c_ 

Fc 


(2.4.1)  and  (2.4.2)  imply  that  limsup  <5q  (x^  ) 

j  -►oo  1 


and  limsuD  Z (x^  )  are  both  finite. 


Suppose  0  is  a  limit  point  of  {Z(x^  )} 

1 


★ 

and  lim  Z (x^  )  =0,  where  {T.  }  is  a  subsequence  of  {T.} 

^If  D 

T  -VOO  J  Ir  ^ 


Jk 


By  proposition  2.7.  and  b  =  — ,  we  have 


lim  u(x^  ) 


T  -voo 


TZ(x*  )  -  (n(x*  )  +  1)  c 

=  lim  c  + - llS - - +  vc_ 


(1  +  Z(x*  ))"^^T.  ^  1  -  1 


Z(x.j,  ) 


:'k 


+  z ( X^  ) u 
^k 


lim  lim  c  t  -  (n^llcZ  ^  vcKli?)"^^ 

n-»-oo  z-^0  Z  [  (1+1^)  -  1) 


-  11 


+  Zu. 


Applying  L' Hospital's  rule  twice,  we  obtain 

2  5. 

lin  TZ  -(n+l)cZ  +  vc  I  (1-kv)  -  l  _  Tv  nc 

Z((1+|)”'^^  -  1]  ^ 


which  implies 


-J4- 


* 

lim  u(x^  ) 

■i  -♦•oj  ^  Ic 


lim 


(n+1) c 
2 


Tv 


n+1 


,  a  contradiction  with  (2.4.2).  Thus, 
(2.4.3)  lim  Z (x*  )  =  Z®  >  0. 

j^-voo 

By  proposition  2.7.  we  have 


QO 


(2.4.4) 


Jk 


)  = 


Z  (Xy  )T  -  (n(x^  )  +  l)c 

jk  ^  ^  ^k 

(1+2  (X*  ))"^^T.  ^  ^-1 


Z  ) 
^k 


Because  (2.4.3)  holds  it  can  be  seen  that  the  first  term 

of  the  right-hand  side  of  (2.3.3)  tends  to  zero  when 
^  c 

thus,  lim  ^^(x^  ^  “  “0  ' 

^k  ^ 


Since  c  +  v— _  +  yZ  >  c  +  v/\ic  +  u  —  is  always  true 

2^  /  V  /yc 

/  V 

for  Z®  >  0,  we  have 


lim  u(x^  ) 
^k  ^k 


*  * 

=  lim  c  +  v  ^^(x^  )  +  yZ(x^  ) 

^k  ^k 


c  „0 

=  c  +  V  ^  +  yZ 


>  c  +  V  /yc  +  y  —  , 

V  ^ 

/  V 
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together  with  (2.4.2)  and  the  definition  of  {x^  },  we  have 


lim  u(x^  ) 

-1  ->oo  -t 


/ —  c 

c  +  v/uc  +  u  — 


/uc 

/  V 


lim  u(x  ) 
1-^“  j 


Q.E.D. 


2.5.  Optimization  by  an  algorithm. 

The  objective  min  g{b,n(b))  in  proposition  2.8.  can  be 
b>0 

standerlized  as  follows. 

min  g (b,  n  (b) ) 
b^O 

bl^_ln(b)  +1]  1  TV  . 

=  C  min  11  +  - --Tf-v-ipi -  +  f-  +  b=T  —  } 

b>0  -1  ^  ^ 


=  C  min  g(b,n(b),~,  s) 
b-0  ^ 


and  by  multiplying  6.’s  by  —  respectively,  (2.3.6) 

X  c 


becomes 


T  V  .  I""’)  "  .  1 

~  c  ^0  ^i_^_b)nfb)+l  _  ^  b 


s^(b]  =  ^  <s^  =  (1+b)^  «Q(b)  - 


for  r=l , 2 , . . . , n (b) . 

—  •  M 

Note  that  q  and  (S^(b)'s  are  determined  by  ^  and 
_  *  * 

b,  and  the  optimal  b  for  g  is  also  optimal  for  g. 

* 

For  the  optimal  b  we  have  following  algorithm. 
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Algorithm  2.1. 


Step  1.  Use  the  bisection  method  to  obtain  the  unique 

zero  of  f  ,,  f  10/®),  lb  , , /®) , 

0  0  0 

[b^  •  t  resnectively . 

If  bj^^  ^  then  n=k  and  go  to  step  2. 

Tv 

step  2.  If  — -  s(b  )  ^  6 ,  go  to  step  3;  otherwise,  let 

^  "o+i 

Tv  “ 

m=max  Ik  ';  —  <  6 ,  n^  <  k  ^  n} ,  go  to  step  4 

Step  3.  If  ^(b  ,T,n.)  •  0  then  use  the  bisection  method 

db  0 

to  obtain  the  unique  zero  ^ 

[c,  b  where  c  is  a  small  positive  number 

such  that  ,n(3)  V  0.  If  <  0, 

let  b  =  b  +1.  Set  m=n_  then  qo  to  step  5. 

"O  "O  ° 


Step  4.  Graph  g(b,n(b))  from  o  to  bm+1,  let  b^  satisfy 

gl^m)  =  min  {g(b,n(b))}.  G(>  to  step  5. 

0<b<b  ,, 

—  —  m+1 


Step  5.  Let  k=m,...,n. 

^  ^*^k+l'*^^  ^  then  use  the  bisection 

'V#  do 

method  to  obtain  the  unique  zero  bj^  of  ^  (‘rk)  in 

then  6,^  =  b|^.  If  ^Ib^.h)  ag<Vl''"’  i  “ 
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0  then  c,(. 

•k 

step  6.  The  optimal  b  satisfies 

*  *  f\. 

q  (b  ,  (n  (b  )  )  =  min 

k=m,  .  .  . ,  ri 


Tl^  stop  6. 


n  (bj^)  )  1  . 


The  finiteness  of  step  1. 

Suppose  lim  b  <  /'c”  ,  by  the  definition  of  b,  and 
k  ^ 

proposition  2.11  we  have 


:k/^ 

/  uv 


Tv  c 
C  MV 


+  1 


0  as  k- 


which  is  false. 


Hence, 


there  is  a  n+1  such  that  b-  , 

n+i 


The  optimization  of  algorithm  2.1.  is  then  followed 
by  the  finiteness  of  stei>  1.,  theorem  2.4,  t>ropo.s  i  t  i  on  2.7 
and  proposition  2.8. 
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